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We study the failure properties of fiber bundles when continuous rupture goes on due 
to the application of external load on the bundles. We take the two extreme models: 
equal load sharing model (democratic fiber bundles) and local load sharing model. The 
strength of the fibers are assumed to be distributed randomly within a finite interval. 
The democratic fiber bundles show a solvable phase transition at a critical stress (load 
per fiber). The dynamic critical behavior is obtained analytically near the critical point 
and the critical exponents are found to be universal. This model also shows elastic-plastic 
like nonlinear deformation behavior when the fiber strength distribution has a lower cut- 
off. We solve analytically the fatigue-failure in a democratic bundle, and the behavior 
qualitatively agrees with the experimental observations. The strength of the local load 
sharing bundles is obtained numerically and compared with the existing results. Finally 
we map the failure phenomena of fiber bundles in terms of magnetic model (Ising model) 
which may resolve the ambiguity of studying the failure properties of fiber bundles in 
higher dimensions. 

Keywords: fiber bundle; critical behavior; universality; fatigue 

1. Introduction 

Fracture processes in heterogeneous media, initiated by external loading has rich 
physical and mathematical aspects^. Since past decades, people are trying to study 
this fracture phenomena through different models as well as through direct experi- 
ments. The knowledge of the strength of a disordered solid and prior knowledge of 
its failure properties are of extreme importance in architectural engineering, textile 
engineering and in construction of any mechanical structure. This knowledge is also 
required to design material microstructures which can be used to construct highly 
reliable components. In the context of geophysics the dynamical aspects of fracture 
process is useful to explain the earthquake phenomena. Some theoretical models 
like lattice models in various dimensions have been proposed to analyse the details 
of the fracture phenomena; among these the fiber bundle model is the earliest and 
the simplest one. 

A loaded bundle of fibers represents the various aspects of fracture process 
through its self-organised dynamics. The fiber bundle model study was initiated 
by Peirce^ in the context of testing the strength of cotton yarns. Since then this 
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model has been studied from various points of view. Fiber bundles are of two classes 
with respect to the time dependence of fiber strength: The 'static' bundles contain 

fibers whose strengths are independent of time, whereas the 'dynamic' bundles are 
assumed to have time dependent elements to capture the creep rupture and fatigue 
behaviors. According to the load sharing rule, fiber bundles are being classified into 
two groups: Equal load-sharing (ELS) bundles or democratic bundles and local load- 
sharing (LLS) bundles. In democratic bundles intact fibers bear the applied load 
equally and in local load-sharing bundles the terminal load of the failed fiber is given 
equally to all the intact neighbors. The classic work of Daniels^ on the strength of 
the static fiber bimdles under equal load sharing (ELS) assumption initiated the 
probabilistic analysis of the model** The distribution of burst avalanches during 
fracture process is a marked feature of the fracture dynamics and can be observed 
in ultrasonic emissions during the fracture process. It helps characterizing different 
physical systems along with the possibility to predict the large avalanches. From a 
nontrivial probabilistic analysis, Hemmer and Hansen^ got power law distribution 
of avalanches for static ELS bundles, whereas the power law exponent observed nu- 
merically for static LLS bundles differs significantly. This observation induces the 
possibility of presenting loaded fiber bundles as earthquake models^. The recent 
mean field estimate of the 'avalanches' in the ELS bundles^'^^, gives a new power 
law"'^'^. The phase transition^'® and dynamic critical behavior of the fracture process 
in such bundles has been established through recursive formulation^'*^'*^'*^ of the 
failure dynamics. The exact solutions*°'*^'*^ of the recursion relations suggest uni- 
versal values of the exponents involved. Attempt has also been made^^ to study the 
ELS ans LLS bundles from a single framework introducing a 'range of interaction' 
parameter which determines the load transfer rule. 

Coleman^^ started working on time dependent bundles under ELS to obtain their 
life time with steady load. Later, subsequent generalization was made by Phoenix^^ 
where both ELS and LLS were considered. Some recent developments*^ show 
fatigue behavior of ELS bundles considering fluctuation in applied load. Also, in- 
troducing noise-induced failure probability of fibers, fatigue behavior is achieved*^ 
in a homogeneous fiber bundle under ELS with steady load. 

In this report, we give a brief summary of the exactly solvable*"'**'*^ static 
fiber bundle models under ELS assumptions. This gives the mean field behavior of 
the failure dynamics, its critical behavior and its universality. An analytic study of 
fatigue-failure in homogeneous bimdlc under ELS*^ is also discussed. Some mod- 
ifications and discussions have been added here to correlate the sequential devel- 
opments. We also show numerically that critical strength of such linear bundles 
under LLS, vanishes in the large chain limit. This is basically a confirmation of 
earlier observations^^'^^. A brief magnetic mapping model of fiber bundles is also 
discussed later. 
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2. Equal load sharing (ELS) bundles 



Fig. 1: The fiber bundle consists initially of N fibers attached in parallel to a fixed plate at the 
top and a movable plate from which a load F is suspended at the bottom. In the ELS model the 
load F is equally shared by all the intact fibers. 

The fiber bundle consists of N fibers, each having identical spring constant K. The 
bundle supports a load F = Na and the breaking threshold (crt/j) ■ of the fibers 
are assumed to be different for different fiber (i). ELS model assumes equal load 
sharing, i.e., the intact fibers share the applied load F = Na equally where a is 
the initial applied stress (load per fiber). The strength of each of the fiber {<Tth)i 
in the bundle is given by the stress value it can bear, and beyond which it fails. 
The strength of the fibers are taken from a randomly distributed normalised density 
p{<7th) within the interval and 1 such that 

/ p{(^th)duth = 1- (1) 
Jo 

The equal load sharing assumption neglects 'local' fluctuations in stress (and its 
redistribution) and renders the model as a mean-field one. 



2.1. Breaking dynamics of the ELS models 

The breaking dynamics starts when an initial stress a (load per fiber) is applied on 
the bundle. The fibers having strength less than a fail instantly. Due to this rupture, 
total number of intact fibers decreases and rest of the (intact) fibers have to bear 
the applied load on the bundle. Hence effective stress on the fibers increases and 
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this compels some more fibers to break. These two sequential operations, the stress 
redistribution and further breaking of fibers continue till an equihbrium is reached, 
where either the surviving fibers are strong enough to bear the apphed load on the 
bundle or ah fibers fail. 

This self organised breaking dynamics can be represented by recursion relations 
in discrete time steps. Let Ui be the fraction of fibers in the initial bundle that 
survive after time step t, where time step indicates the number of occurrence of 
stress redistribution. Then the redistributed load per fiber after t time step becomes 



and after t+1 time steps the surviving fraction of fiber is 

[/t+i =1-P(at); (3) 
where P{ot) is the cumulative distribution of the corresponding density p{oth)'- 

P{at) = [ * p{ath)dath. (4) 
Jo 



Now using Eq. (2) and Eq. (3) we can write the recursion relations which show how 
at and Ut evolve in discrete time: 



and 



Ut+i = l-P{a/Uty,Uo = l. (6) 



The recursion relations (5) and (6) represent the basic dynamics of failure in 
equal load sharing models. At the equilibrium or steady state Ut+i = Ut = U* and 
at+i = at = a* . This is a fixed point of the recursive dynamics. Eq. (5) and Eq. 
(6) can be solved at the fixed point for some particular form of p{ath) and these 
solutions near U* (or a*) give the detail features of the failure dynamics of the 
bundle. 
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2.2. Phase transition and critical behavior for uniform 
distribution of fiber strength 




Fig. 2: The simplest model considered here assumes uniform density p{(Tth) oi the fiber strength 
distribution up to a cutoff strength. 

We choose the uniform density of fiber strength distribution to solve the recursive 
failure dynamics of democratic bundle. Here, the cumulative distribution becomes 

P{(Jt) = / p{<Jth)dath = I dath = crt- (7) 
Jo Jo 

Therefore Ut follows a simple recursion relation (following Eq. (6)) 

C^m = l-^. (8) 

At the equihbrium state (C/t+i = Ut = C/*), the above relation takes a quadratic 
form of U* : 

U*^ -U* + a = 0. (9) 

The solution is 

[/*(<7) = i±(ae-c7)V2;^^=l. (10) 

Here ctc is the critical value of initial applied stress beyond which the bundle fails 

completely. The solution with (+) sign is the stable one, whereas the one with (— ) 
sign gives unstable solution^-^'^^. The quantity U*{a) must be real valued as it has 
a physical meaning: it is the fraction of the original bundle that remains intact 
under a fixed applied stress a when the applied stress lies in the range < a < ac- 
Clearly, U*{(Tc) = 1/2 (putting cr = CTc in Eq. 10). Therefore the stable solution can 
be written as 

U*{a) = C/*(ae) + (ae - a)^/^; = ]. (11) 
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For a > ac we can not get a real-valued fixed point as the dynamics never stops 
until Ut = when the bundle breaks completely. It may be noted that the quantity 
U*{a) — U*{(Jc) behaves like an order parameter that determines a transition from 
a state of partial failure (fj < ctc) to a state of total failure {a > (Jcy^'^"^- 

= U*{a)-U*{a,) = {a,-af;(i=l. (12) 



(a) At cr < Ctc 

To study the dynamics away from criticality (ct — > CTc from below), we replace the 
recursion relation (8) by a differential equation 

_dU_lP-U + a 

dt~ u ■ ^ ^> 

Close to the fixed point we write Ut{a) = U*{a) + e (where e ^ 0). This, following 
Eq. (10), gives^O'" 

e = Ut{a) - U*{a) « exp(-VT), (14) 
where t = ^ [|(ctc — ct)~^/^ + l] . Near the critical point we can write 

T oc ((Tc - o-)"";a = ^- (15) 

Therefore the relaxation time diverges following a power-law as ct ^ CTc from 
belowi°'i^ 

One can also consider the breakdown susceptibility %, defined as the change of 
U*{a) due to an infinitesimal increment of the applied stress o-^O'ii'i^ 



dU*ia) 



da 



i(CT,-CT)-^;7=^ (16) 



from equation (10). Hence the susceptibility diverges as the applied stress a ap- 
proaches the critical value ac = j- Such a divergence in x had already been observed 
in the numerical studies'"'^. 

(b) At (7 = a"c 

At the critical point (ct = CTc), we observe a dynamic critical behavior in the relax- 
ation of the failure process to the fixed point. From the recursion relation (8) it can 
be shown that decay of the fraction Ut{<Jc) of unbroken fibers that remain intact at 
time t follows a simple power-law^ 

+ (17) 

starting from Uq = 1. For large t [t ^ oc), this reduces to Ut — 1/2 oc t^^; 5 = 1; a 
strict power law which is a robust characterization of the critical state. 
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2.3. Universality class of the model 

To check the universality class of the model we have taken two other types of fiber 
strength distributions: linearly increasing density distribution and linearly decreas- 
ing density distribution within the limit and 1. We solve the recursion equations 
(5) and (6) in these two cases. We show that while ctc changes with different strength 
distributions, the critical behavior remains unchanged: q; = 1/2 = /3 = 7, 5=1. 

(a) Linearly increasing density of fiber strength 




Fig. 3: The linearly increasing density p{<Tth) of the fiber strength distribution up to a 

cutoff strength. 

Here, the cumulative distribution becomes 

P{<^t) = / p{(^th)d(7th = 2 / (Jthdoth = <y1- (18) 
Jo Jo 

Therefore Ut follows a recursion relation (following Eq. (6)) 

U,,. = l-[^y. (19) 

At the fixed point {Ut+i = Ut = U*), the above recursion relation can be represented 
by a cubic equation oi U* 

{U*f-{U*f + a^ = 0. (20) 

Solving the above equation we get the value of critical stress cTc = a/4/27"'"^ 
which is the strength of the bundle for the above fiber strength distribution. Here, 
the order parameter can be defined as U*{a) — U*{ac) and this goes as 

0^{ae-af;P=^. (21) 
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The susceptibility diverges as the critical point is approached from below: 



X 



dU*ia) 



da 



oc (cTc -a) = 



1 



We can also show that for any a < ac 

Ut{a)-U*{a)^eM-t/T), 

with 



and at a = ac 



T oc {ac — a) " ; a = 



Ui - U*{ac) oc t-^:6= 1. 



(b) Linearly decreasing density of fiber strength 




(22) 

(23) 
(24) 

(25) 



Fig. 4: The linearly decreasing density p(<7th) of the fiber strength distribution up to a cutofT 

strength. 



In this case, the cumulative distribution becomes 

/■ft pert 

P{(^t) = / P{(^th)dath = 2 / (1 - ath)dath = 2at - a^ 
Jo Jo 

and Ut follows a recursion relation (following Eq. (6)) 

At the fixed point (C/t+i = Ut = U*), the above recursion relation can be represented 
by a cubic equation of U* 

{U*f-{U*f+2aU*-a'^ = 0. (28) 



(26) 



(27) 
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Solution of the above equation suggests the value of critical stress CTc = 4/27^^ 
which is the strength of the bundle for the above fiber strength distribution. Also, 
the order parameter goes as 



O = [U*ia) - U*{a,)] oc {a^ - af;f3= ]- 



(29) 



and the susceptibility diverges with the similar power law as in Eqs. (16 and 22) 
when the critical point is approached from below: 



X 



Here also for any a < Uc 



where 



and sX a = Ur 



dU*{a) 



da 



oc ((7c -a) : 



Ut{a)-U*{a)^eM-tM, 



Ut-U*{ac)oc.t-^\5 = l. 



(30) 

(31) 
(32) 
(33) 



Thus the democratic fiber bundles (for different fiber strength distributions) 

show phase transition with a well defined order parameter which shows similar power 
law variation on the way the critical point is approached. The susceptibility and 
relaxation time also diverge with same power exponent for all the cases. Therefore, 
failure of democratic fiber bundles belong to a universality class characterized by 
the universal values of the associated exponents (a, (3, 7 and 5). 



2.4. Nonlinear stress-strain relation 



One can now consider a slightly modified strength distribution of the democratic 
fiber bundle, showing nonlinear deformation characteristics^'^'^^. For this, we con- 
sider an uniform density distribution of fiber strength, having a lower cutoff. Until 
failure of any of the fibers (due to this lower cutoff), the bundle shows linear elastic 
behavior. As soon as the fibers start failing, the stress-strain relationship becomes 
nonlinear. The dynamic critical behavior remains essentially the same and the static 
(fixed point) behavior shows elastic-plastic like deformation before rupture of the 
bundle. 

Here the fibers are elastic in nature having identical force constant k and the 
random fiber strengths distributed uniformly in the interval [ctl, 1] with ctl > 0; 
the normaUsed distribution of the threshold stress of the fibers thus has the form 
(see Fig. 5): 

p{(Tth) = < 1 ^ ^ ^ < 1 • 34 
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Fig. 5; The fiber breaking strengtii distribution p{<Jth) considered for studying elastic-plastic type 
nonlinear deformation behavior of the ELS model. 

For an applied stress u < gl none of the fibers break, though they are elongated 
by an amount 5 = a/n. The dynamics of breaking starts when applied stress a 
becomes greater than ctl. Now, for a > ul the fraction of unbroken fibers follows a 
recursion relation (for p{(Jth) as in Fig. 5): 



Ut+i = 1 - 



NUt 

which has stable fixed points: 
1 



2(1 -ai) 



CTL 



1 + 



1-<JL 



1-^ 

CTc 



1-^ 

Ut 



1/2 



4(1 -(Tl)' 



(35) 



(36) 



The model now has a critical point CTc = 1/[4(1 — ctl)] beyond which total failure 
of the bundle takes place. The above equation also requires that (Jl < 1/2 (to keep 
the fraction U* < 1). As one can easily see, the dynamics of Ut for a < ac and also 
at (7 = (Jc remains the same as discussed in the earlier section. At each fixed point 
there will be an equilibrium elongation ^(fT) and a corresponding stress S = U*KS{a) 
develops in the system (bundle). This 5 (a) can be easily expressed in terms of U*{a). 
This requires the evaluation of a*, the internal stress per fiber developed at the fixed 
point, corresponding to the initial (external) stress a (= F/N) per fiber applied on 
the bundle when all the fibers were intact. From the first part of Eq. (36), one then 
gets (for a > ctl) 



(1 - Cl) 1- <7L 



Consequently, 



KS{a) = a* = 1 - (1 - aL)U*{a). 



(38) 
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It may be noted that the internal stress a* (= ac/U*{ac) from (2)) is universally 
equal to 1/2 (independent of cri;,;from (36)) at the failure point a = CTc of the bundle. 
This finally gives the stress-strain relation for the ELS model : 




K 2k 



Fig. 6: Schematic stress (S)-strain (<5) curve of tlic bundle (shown by the solid line), following Eq. 
(39), with the fiber strength distribution (34) (as shown in Fig. 5). 

This stress-strain relation is schematically shown in Fig. 6, where the initial linear 
region has slope k (the force constant of each fiber). This Hooke's region for stress S 
continues up to the strain value 5 = (tl/k. until which no fibers break {U*{a) = 1). 
After this, nonlincarity appears due to the failure^ of a few of the fibers and the 
consequent decrease of U*{a) (from unity). It finally drops to zero discontinuously 
by an amount a*U*{ac) = 1/[4(1 — ctl)] = <7c at the breaking point cr = CTc or 
5 ^ (j*/k = 1/2k for the bundle. This indicates that the stress drop at the final 
failure point of the bundle is related to the extent {ul) of the linear region of the 
stress-strain curve of the same bundle. 
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2.5. Fatigue-failure in a homogeneous fiber bundle 



Fatigue in fiber bundle model was first studied by Coleman in 1958^^. Thermally 
activated failures of fiber have recently been considered and approximate fatigue 
behavior has been studied^^. We consider here a very simple fiber bundle model 
with noise-induced activated failure, for which the dynamics can be analytically 
solved. 

Let us consider a homogeneous bundle of N fibers under load F{= Na), each 
having identical failure strength ath which is the strength (ctc) of the bundle also. 
Without any noise (T = 0), the model is trivial: the bundle does not fail (failure 
time T is infinity) for stress a < Uc and it fails immediately (r = 0) for a > Gc- We 
now assume that each such fiber has a finite probability p{cF, T) of failure at any 
stress <j induced by a non-zero noise T: 



;;7exp 



1, 



< cr < Ctc 
CT > Ctc 



(40) 



As one can see, each fiber now has got a non-vanishing probability p{a, T) to fail 
under a stress a < a,, at any non-zero noise parameter T. p{<T, T) increases as T 
increases and for ct > Uc, p{<J,T) = 1. Unlike at T = 0, the bundle therefore fails at 
a < ac after a finite time r. Here we assume each fiber to have unique threshold, 
while their breaking probability at any a (< CTc) is due to noise-activated hopping 
over the barrier height (ctc — a). This differs from the earlier model studies^^'^^ 
where the load distribution is noise induced. 

Failure time 

At T ^ and under any stress ct (< CTc), some fibers fail due to noise and the load 
gets shared among the surviving fibers, which in turn enhances their stress value, 
inducing further failure. Denoting the fraction of fibers to the initial bundle that 
remains intact at time t by Ut, a discrete time recursion relation^^ can be written 
as 



Ut+i = Ut 



1-p 



(41) 



where a/Ut = Fj {NUt) is the redistributed load per fiber among the NUt surviving 
fibers at time t. In the continuum limit, we can write the above recursion relation 
in a differential form 

dU 

= — exp - 

T 



dt 



CT 

= — exp 

CTc 



-ir-u-i) 

r ^ CT / 



(42) 



The failure time r is defined as t = t when Ut = 0. Integrating Eq. (42) within 
proper limits we get 

r-l 



= / dt=—exp( — -L\ [ exp i f — 

Jo 0- \ tJ Jo Lr ^ 



u 



dU 



(43) 



I 



I 
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or 



17 



(44) 

for a < (Tc- For a > (Tc, starting from Ut = I aX t = Q, one gets Ut+\ = from Eq. 
(41), giving T = 0. 



For small T and as a 



(ae/(7-l)/T 



. This failure time r 



Tc, T ~ Texp 

therefore approaches infinity as T — > 0. For cr < CTc, one gets finite failure time r 
which decreases exponentially as a approaches (7^ oi' T increases and r = for 
a > (Tc- This last feature is absent in the earlier formulations^®. However, all these 
features are very desirable and are in qualitative agreement with the recent ex- 
perimental observations^^. Our numerical study confirms the above analytic results 
[obtained using the continuum version of the recursion relation (41) (see Fig. 7) well. 



(a) 



3 OS 




200 400 600 too 1000 1200 1400 1600 
t 




0.12 0.14 0.16 0.18 



Fig. 7. The simulation results showing variation of average failure time t against (a) stress a and 

(b) against noise T, for a bundle containing A'' = 10^ fibers. The theoretical results are shown by 
dotted and dashed lines [from Eq. (44)]. The insets show the simulation results for the variation 
of the fraction U of unbroken fibers with time t for different T values [1.2 (cross) and 1.0 (plus)] 
in (a) and a values [0.15 (cross) and 0.12 (plus)] in (b). The dotted and dashed lines represent the 
theoretical results [Eq. (42)]. 



3. Strength of the local load shciring (LLS) fiber bundles 



The strength of a bundle of fibers plays important role in the failure dynamics 
of the bundle when external load is applied. By a probabilistic approach Daniel'^ 
pioneered the study of finding the strength distribution of a bundle in terms of 
strength of the constituent fibers. Daniel assumed equal sharing of applied load 
(ELS). Later, this work has been expanded in the context of material science^. 
This type of model shows (both analytically and numerically) existence of a critical 
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strength (non zero Uc) of the bundle®'^"'^^'^^ beyond which it collapses instantly. The 
other extreme model, i,e., the local load sharing (LLS) model has been proved to be 
difficult to tackle analytically^^. However some approximate asymptotic mcthods^-'^ 
have been developed to tackle the problem in one dimension. Recently Pacheco 
et. al.^^ introduced and solved the one sided load transfer model analytically. It is 
basically a simplification of the conventional (both sided) LLS model. Considering 
WcibuU distribution of fiber strength, they obtained the system size dependence of 
the strength (ctc) of the bundle as 

lim — = const + logj N. (45) 

JV— >oo fJc 

The subscript s indicates the shape factor or WeibuU index of the fiber strength 
distribution. This clearly shows cTc — > as A/' — > oo. Smith^^ conjectured a similar 
logarithmic dependence for LLS bundles from numerical results. 

Here, we have simulated the above two types of LLS models: the one sided load 
transfer model and the conventional both sided load transfer model considering 
uniform distribution (random) of fiber strengths [Fig. 8] . We can not use the strictly 
uniform distrimition (as in Fig. 2), to avoid the increasing sequence of fiber strength 
arrangement which is fatal in case of LLS bundles though it does not matter in case 
of ELS bundles. 

1.5 r 



1.5 



Fig. 8: The uniform (random) fiber strength distribution p{(7th) considered to study the LLS 

models. 
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We observe the 1/ log dependence of the bundle's strength (ctc) for both the cases 
[Fig. 9] which confirms the non existence of any critical strength (non zero CTc) of 
the bundle in one dimension. 

16 
15 
14 
13 
12 

--ill 
»— I 

10 

9 
8 
7 
6 

100 1000 10000 100000 

N 

Fig. 9: The system size dependence of strength of the bundle ((7c) in LLS models. The results are 
obtained after averaging over 5 x lO^configurations. 

4. Magnetic mapping of the fiber bundle models 

We can map the fiber bundle problem in terms of a magnetic model or specifically 
an Ising model. The state of the fibers represent different states of the spins. We 
introduce a variable Ui{t) to indicate the status of the i-th fiber at the time step t 
such that Ui (t) — 1 for intact fibers and for broken fibers. The corresponding spin 
variable Si{t) = 2ui{t) — 1 can then take two values ±1, which can represent the up 
and down states of the Ising spin at i-th site. The strength threshold {ath)i of the 
i-th fiber represents a random field at the i-th site (fiber) favoring Si = +1 state (or 
Ui = 1) and the stress ai{t) on the i-th fiber represents a competing field at each 
site favoring the state Si = —1 (or m =0). The value of Si{t) or Ui{t) depends on 
the strength of the resultant field hi{t): 

Si{t)=sgn[hi{t)]; (46) 

hi{t) = {ath)i - cTi{t);Ui{t) = ^ [si{t) + 1] . (47) 
The only restriction on the competing fields is 
^ 1 

Y,Uiit)(Ti{t) = - [s,(i) + 1] ai{t) = F; (48) 
i=i ^ 
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F is the applied force (load) on the bundle. ai{t), and hence hi{t), evolves with time 
t following the above relation. This is a general restriction for both, equal and local 
load sharing cases. The average magnetization at time t can be defined as 

1 ^ 

m{t) = j;^J2'^(ty (49) 

i=l 

This non-interacting spin model is somewhat special as the zcro-tcmpcrature dy- 
namics of the spins here is determined by the strength of the effective fcild hi at 
that site (46), where one part {(Ti{t)) of hi evolves with time following (48). The 
stable manetization m (= 1 at cr = 0) decreases continuously with a until it reaches 
a value nic = 2U* — 1, where 

1 ^ 

i7* = -^w.(t^oo), (50) 

at a = <7c, beyond which m drops to —1 discontimiously. 
For equal load sharing case, the above equation becomes 

N 

a{t)J2ui{t) = F = aN; (51) 

i=l 

a being the initial applied stress (F/N). For uniform fiber strength distribution we 
have U*{ac) — 1/2 (from (10)). Here, rric = and m decreases continuously from 
unity (at (7 = 0) to m = rric = at (7 = (Tc and then abruptly crosses over to 
m = —1 at cr > CTc. 

5. Discussions 

The inherent mean-field nature of the ELS models enables to construct recursion 
relations (Eqs. 5 and 6) which captures essentially all the intriguing features of the 
failure dynamics. Though we have identified O = U*(a) — U*{ac) oc ((Jc — o")^ as 
the order parameter (with exponent /3 = 1/2) for the continuous transition in the 
ELS models, unlike in the conventional phase transitions it does not have a real- 
valued existence for a > ac- The 'type' of phase transition in ELS models is still 
a controversial issue. Earlier^ it was suggested to be a first order phase transition, 
because the the surviving fraction of fibers has a discontinuity at the breakdown 
point of the bundles. However, as the susceptibility shows divergence (x oc (uc — 
'jI = 1/2) at the breakdown point, the transition has been later identified 
to bo of second order*'*''^^'^^. The dynamic critical behavior of the ELS models 
and the universality of the exponent values are straightforward. Here, divergence of 
relaxation time (r) at the critical point (r oc (ctc — cr)~"; a = 1/2) indicates 'critical 
slowing' of the dynamics which is characteristic of conventional critical phenomena. 
At the critical point, one observes power law decay of the surviving fraction in time 
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(t^t(o'c) oc t~^;d = 1). We demonstrated the universality of the failure behavior 
near a = a^, for three different distributions: uniform (Fig. 2), linearly increasing 
(Fig. 3) and linearly decreasing (Fig. 4) distributions of fiber strength. The critical 
strengths of the bundles differ in each case: CTc = 1/4, \jA/27 and 4/27 respectively 
for these three distributions. However, the critical behavior of the order parameter 
O, susceptibility x, relaxation time r and of the time decay at (Tc, as given by the 
exponents /3, 7, a and 5 remain unchanged: a = 1/2 = /? = 7 and (5 = 1 for all three 
distributions. 

The ELS model also shows realistic nonlinear deformation behavior with a 

shifted (by (7^, away from the origin) uniform distribution of fiber strengths. The 
stress-strain curve for the model clearly shows three different regions: clastic or 
linear part (Hooke's region) when none of the fibers break {U*{a) = 1), plastic 
or nonlinear part due to the successive failure of the fibers {U*{a) < 1) and then 
finally the stress drops suddenly (due to the discontinuous drop in the fraction of 
surviving fibers from U*{ac) to zero) at the failure point cTc = 1/[4(1 — ctl)]- 

The fatigue study in a homogeneous fiber bundle suggests if the each fiber has a 
finite probability of failure (due to noise T as in (40)) below its normal strength, then 
the failure time of the bundle decreases exponentially (t ~ Texp ((Tc/(t — 1) /T ) 
as cr approaches cTc from below and r ~ for a > cTc. These features agree well with 
the experimental observations in disordered solids. 

The LLS bundles show 'zero' critical strength as the bundle size goes to infinity 
in one dimension. It is not clear at this stage if, in higher dimensions, LLS bundles 
are going to have non-zero critical strength^ In any case, the associated dynamics 
of failure of these higher dimensional bundles with variable range load transfer^^'^^ 
should be interesting. 

We believe, the elegance and simplicity of the model, its common-sense appeal, 
the exact solubility of its critical behavior in the mean field (ELS) limit, its demon- 
strated universality, etc, would promote the model eventually to a level competing 
with the Ising model of magnetic critical behavior. 
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